AD-A061  462  CALIFORNIA  UNIV  SAN  DEI&O  LA  JOLLA  DEPT  OF  ELECTRICA— ETC  F/6  3/2 

THE  LATITUOINAL  STRUCTURE  OF  SOLAR  HIND  STREAMS  FROM  RADIO  SCIN--ETC(U) 
SEP  79  S  ANANTHAKRISHNAN#  H  A  COLES  F1962S-77-C-0161 

UNCLASSIFIED  SCIENTIFIC-2  AFGL-TR-79-0234  NL 


l 

1 

i 

'4- 

END 

4^80 

MICROCOPY  RESOLUTION  TEST  CHtfU 
national  bureau  or  standards-  1%3-jl 


FILE  COPYi  AD  A  0  8 1  4  8  '<& 


AFGL-TR-79-0234 


THE  LATITUDINAL  STRUCTURE  OF 
SOLAR  WIND  STREAMS  FROM 
RADIO  SCINTILLATION  OBSERVATIONS 

S.  Ananthakrlshnan 
W.A.  Coles 
J.J.  Kaufman 
B.J.  Rlckett 


Department  of  Electrical  Engineering 
and  Computer  Sciences 
University  of  California,  San  Diego 
La  Jolla,  California  92093 


<3 y 


Scientific  Report  No.  2 


September  1979 


Approved  for  public  release;  distribution  unlimited 


AIR  FORCE  GEOPHYSICS  LABORATORY 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
HANSCOM  AFB,  MASSACHUSETTS  01731 


s 


DT1C 

.electei 


MAR  6  I960 


80 


4  004 


REPORT  DOCUMENTATION  PAGE 


J’HE  LATITUDINAL  STRUCTURE  OF^OLAR^ 
WIND  STREAMS  FROM  jlADIO  SCINTILLATION 

Observations  '  — - - ^ 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


6.  PERFORMING  ORG.  REPORT  NUMBER 


— -J — - — ^ ^ - J 

S./Ananthakrishnan  B.J./Rickett  / 

W.Ay/Coles  — ~~~ - 1 

J0  J.  /Kaufman  j 


9.  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS  / 

Dept  of  Electrical  Engineering  and  Computer  Sci< 
University  of  California,  San  Diego 
La  Jolla,  California  92093 


II.  CONTROLLING  OFFICE  NAMe  ANO  ADORE** 

Air  Force  Geophysics  Laboratory  *j 

Hanscom  AFB,  Massachusetts  01731  JL*' 

Monitor/ Richard  C.  Altrock/ PHS 


F 1 96  2  8  -  77  -C  -ffl  6 1 


U.  MONITORING  AGENCY  NAME  A  ADORES S(lt  different  frojfl  Controlling  Office)  IS  SECURITY  CLASS,  (of  thle  report) 


Unclassified 


IS*.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  (hit  Report) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  abatrmct  entered  in  Block  20,  ft  different  from  Report) 


19.  KEY  WORDS  (Contlmi*  on  revree  tide  if  neceaamry  and  identify  by  block  number) 


20.  ABSTRACT  (Conllnu*  on  reverae  aide  if  neceaamry  end  Identity  by  block  number) 


We  present  a  comparison  of  spacecraft  and  IPS  measurements  of  the  solar 
wind  velocity  for  1973-75,  when  the  large  scale  velocity  structure  was  very 
stable.  From  this  comparison  we  infer  a  prototype  distribution  across  a 


DD  ,j2nM73  1473  coition  of  i  nov  65  is  obsolctc  Unclassified  L/ f  }  (qc^/ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  D«U  Entered) 


SECURITY  cl ASS1FICATION  OF  THIS  PAOgfHTifi  Dmt*  -  - 

high  speed  stream  for  the  electron  density  fluctuations,  Anfi,  in  the  spatial 

frequency  range  10”^  -  10”^  km”^  .  The  "microturbulence"  is  more  closely 

related  to  the  velocity  structure  in  a  stream  than  to  the  bulk  density; 

enhancements  in  Ang  occur  wherever  high  speed  wind  overtakes  low  speed  wind. 

The  distribution  of  An  is  closely  matched  by  the  distribution  of  magnetic 

e 

field  fluctuations  ~I0"6  km"^  . 


I 

! 

I 

I 

i 

\ 


Unclassified 


SECURITY  cl  ASSIFtCATlON  OF  THIS  PAOEfWh«n  Enf*r*d> 


Preface 


Since  1972  the  UCSD  74-MHz  radio  observatory  has  continuously  monitored  the 
solar  wind  over  a  wide  range  in  heliographic  latitude  using  the  method  of  inter¬ 
planetary  scintillations  (IPS).  A  primary  goal  has  been  to  study  the  connection 
between  high  speed  streams,  which  cause  geomagnetic  disturbances,  and  the  structure 
of  the  lower  corona.  In  this  report  we  present  in  preprint  form  two  papers  which 
examine  certain  aspects  of  this  investigation.  The  first,  "Microturbulence  in 
Solar  Wind  Streams",  discusses  the  interpretation  of  the  velocity  measured  by  IPS, 
which  is  a  path-smoothed  version  of  the  point  solar  wind  velocity.  This  work 
provides  a  prototype  line-of-sight  smoothing  function  which  can  be  used  to  infer 
point  velocities  from  IPS  data  in  regions  where  spacecraft  measurements  are  not 
available.  The  second  paper,  "Solar  Cycle  Changes  in  the  High  Latitude  Solar  Wind", 
reports  the  discovery  of  a  systematic  variation  with  the  sunspot  cycle  of  the  average 
high  latitude  velocity.  This  variation  is  explained  by  relating  it  to  changes  in  the 
gross  magnetic  field  structure  in  the  lower  corona. 
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"Microturbulence  in  Solar  Wind  Streams" 

S.  Ananthakrishnan,  W.A.  Coles  and  J.J.  Kaufman 


Abstract 


j  comparison  of  spacecraft  and  IPS  measurements  of  the  solar 
—  juia  fcacifiOta.^  /v 


wind  velocity  for  1973-75^  when  the  large  scale  velocity  structure  was  very 
stable.  From  this  comparison  we  infer  a  prototype  distribution  across  a  Jslfe  w 
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enhancements  in(AnN occur  wherever  high  speed  wind  overtakes  low  speed  wind. 
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The  solar  wind  velocity  estimated  from  observations  of  interplanetary 
scintillations  (IPS)  can  be  compared  with  spacecraft  measurements  during 
1973-75  when  the  large  scale  structure  was  very  stable.  Coles,  Harmon, 

Lazarus,  and  Sullivan  (1978)  made  such  a  comparison  using  74  MHz  data  from 
the  ecliptic  source  3C144  both  to  "calibrate"  the  IPS  data  of  1973  and  to  infer 
the  distribution  of  the  electron  density  fluctuations  which  cause  IPS.  In  this 
paper  we  have  estimated  the  distribution  of  this  "micro turbulence"  more  precisely 
using  a  refined  technique,  and  we  have  extended  the  comparison  through  1975. 

The  IPS  observed  by  the  UCSD  74  MHz  observatory  is  caused  by  micro- 

-4  -2  -1 

turbulence  of  spatial  frequency  10  to  10  km  .  The  level  of  the  flue- 
tuation  spectrum  in  this  range,  (Ang)  ,  appears  as  a  weighting  function  in  the 
path  integral  for  an  IPS  observation.  Thus  to  "invert"  this  integral  -  that 
is,  to  recover  the  point  measurement  -  one  must  have  accurate  knowledge  of 
the  distribution  of  (Ang)  .  Unfortunately  Ang  is  not  regularly  observed  by 
spacecraft  because  the  frequency  range  is  too  high,  so  we  must  use  model -fitting 
techniques.  The  IPS  observations  are  reduced  to  three  parameters:  the  scin¬ 
tillation  index,  peak  velocity,  and  midpoint  velocity.  The  index  is  the 
normalized  rms  intensity  variation.  It  is  proportional  to  Ang  in  weak 
scattering,  but  inversely  related  to  Ang  in  strong  scattering.  The  midpoint 
velocity  gives  a  weighted  mean  of  the  velocity  component  transverse  to  the 
line  of  sight,  whereas  the  peak  velocity  selects  the  higher  velocities  along 
this  path.  These  parameters  can  be  calculated  in  weak  scattering  (for 
elongations  greater  than  40°  at  74  MHz)  if  the  velocity  and  Anfi  are  known. 

The  details  of  this  calculation  are  given  by  Coles,  Harmon,  Lazarus,  and 
Sullivan  (1978).  In  strong  scattering  the  velocities  calculated  are  no  longer 
exact  (although  they  are  still  useful),  but  the  index  must  be  discarded. 
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In  the  earlier  work,  the  spacecraft  velocities  and  a  model  distribution 

of  Ang  were  used  to  predict  the  IPS  parameters.  The  predictions  were  compared 

with  the  observations  and  inferences  about  the  model  An  were  drawn.  A  model 

e 

for  which  Ang  “  np  (the  bulk  proton  density)  was  tested.  This  correctly 
predicted  the  increases  of  index  which  are  observed  before  high  speed  streams 
(Houminer  and  Hewish,  1974),  but  the  velocity  agreement  was  poor.  Coles,  Harmon, 
Lazarus,  and  Sullivan  (1978)  found  that  this  model  overestimated  Ang  in  the 
compression  region  and  underestimated  An„  in  the  interior  of  a  stream.  Other 
simple  models  were  tried,  but  none  gave  a  uniformly  good  match.  We  have  confirmed 
these  conclusions  by  extending  the  earlier  analysis  through  1975  (see  also 
Harmon,  1975). 

Methods 

The  distribution  of  microturbulence  can  be  estimated  using  an  iterative 
process  provided,  of  course,  that  it  converges  to  a  unique  solution.  Here  the 
IPS  parameters  are  predicted,  as  before,  and  the  model  An0  is  perturbed  until 
a  "best  fit"  to  the  observed  parameters  is  obtained.  We  have  used  an  inter¬ 
active  program  in  which  an  operator  judges  the  goodness-of-fit  from  a  graphics 
display  and  specifies  the  next  perturbation  on  this  display.  This  was 
necessary  to  speed  convergence  and  to  permit  certain  features  of  the  data,  such 
as  the  rising  edge  of  a  high  speed  stream,  to  be  emphasized.  It  also  permits 
us  to  ignore  doubtful  data  points,  but  makes  the  process  somewhat  subjective. 

We  applied  this  process  to  data  of  1973-75  using  two  fitting  criteria.  First 
we  matched  only  the  peak  and  midpoint  velocities,  ignoring  the  indices.  This 
process  converged,  but  occasionally  two  quite  different  solutions  were  possible. 

We  then  matched  all  three  parameters  simultaneously.  This  additional  constraint 
eliminated  the  ambiguities  and  improved  the  convergence.  It  was  found  that 


introducing  the  index  did  not  degrade  the  agreement  between  predicted  and 
observed  velocities.  Thus  the  index  is  indeed  acting  as  a  constraint,  and 
it  is  consistent  with  the  velocity  data.  We  used  the  index  as  secondary 
data  because  its  sensitivity  to  interference  is  greater  and  it  cannot  be  used 
in  strong  scattering. 

The  prediction  calculation  is  done  exactly  as  described  by  Coles,  Hannon, 

Lazarus,  and  Sullivan  (1978).  It  includes  a  power-law  density  spectrum  with 

an  exponent  of  -2.7  as  described  by  Coles  and  Harmon  (1978);  it  includes  the 
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radio  source  structure;  it  uses  a  solar  distance  dependence  Ang  «  r 
(Armstrong  and  Coles,  1978);  and  it  uses  a  simple  spiral  model  to  map  space¬ 
craft  data  to  the  line  of  sight.  As  before,  no  allowance  is  made  for  stream 
evolution  (with  distance  or  time),  and  heliographic  latitude  differences  of 
up  to  7°  are  ignored.  The  spacecraft  data  used  in  the  calculations  has 
6-hour  resolution,  and  the  iterative  process  is  applied  to  overlapping  30-day 
blocks  of  data.  The  initial  guess  of  A ng  «  np  is  perturbed  first  to  match 
the  velocity,  then  to  match  the  index  without  degrading  the  velocity  match. 

The  final  result  is  an  estimate  of  the  microturbulence,  Ane>  for  1973-75, 
which  can  be  compared  directly  with  spacecraft  plasma  and  magnetic  field 
data. 

Results 

The  spacecraft  velocity  and  density;  the  IPS  parameters;  and  the  Ang 
estimates  for  1973-75  are  shown  in  Figures  1-3.  Only  the  midpoint  velocity 
is  shown,  as  it  is  representative  of  the  velocity  match  and  is  more 
sensitive  to  Ang  than  is  the  peak  velocity.  Spacecraft  (IMP  6,  7  and  8) 
data  gaps,  indicated  by  x's,  were  filled  by  linear  interpolation.  The 
measured  IPS  data,  indicated  by  x's,  are  overplotted  on  the  model  calculation, 
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and  generally  good  agreement  can  be  seen.  For  comparison,  results  for  the 
"initial  guess"  model,  Ang  «  n^,  are  also  shown.  During  this  entire  period 
the  persistence  of  the  overall  velocity  structure  is  strong,  especially  during 
1974  when  two  streams  remain  stable  for  at  least  7  rotations.  However,  there 
are  significant  changes  in  the  streams,  particularly  in  n^  and  An0  at  their 
leading  edges,  even  when  the  velocity  structure  repeats  well. 

A  major  cause  of  error  is  the  "noise"  on  the  IPS  observation.  There  is 
a  random  component  simply  due  to  receiver  noise  and  systematic  component 
caused  by  interference.  These  are  most  important  when  the  scintillations  are 
weak  -  that  is,  when  the  solar  elongation  is  large  (greater  than  120°)  or  small 
(less  than  20°).  Interference  always  increases  the  apparent  scintillation 
index  but  will  only  affect  the  velocities  if  it  is  correlated  over  the  100  km 
spacing  between  antennas.  Solar  radio  bursts  cause  correlated  interference 
with  zero  time  lag,  and  thus  the  velocity  tends  to  be  overestimated;  otherwise 
there  is  no  significant  bias  to  the  velocity,  although  the  errors  do  increase 
as  the  scintillation  decreases.  Typical  estimation  errors  are  10%  for  the 
velocity  and  30%  for  the  index. 

There  are  periods  in  Figures  1-3  where  a  satisfactory  match  could  not  be 
obtained  for  an  Ang  model.  For  example,  near  days  210  and  235  of  1974  the 
velocity  and  index  could  be  matched  separately  but  not  simultaneously.  Here 
a  compromise  is  shown,  but  it  is  clear  that  one  of  our  assumDtions  has 
broken  down.  Occasionally  the  velocity  cannot  be  matched  by  any  model.  For 
example  near  day  80  of  1974  the  IPS  velocities  are  lower  than  any  calculation, 
whereas  near  day  220  of  1974  the  reverse  is  true.  There  are  several  cases, 
such  as  days  113  and  238  of  1973,  where  rapid  changes  occur  in  both  the  space¬ 
craft  and  IPS  velocities  but  do  not  appear  in  the  predictions  (which  have  a 


smoothing  effect).  These  discrepancies  can  be  caused  by  temporal  variations, 
stream  evolution,  or  latitude  differences,  but  they  do  not  affect  a  significant 
fraction  of  the  data. 

The  "best  fit"  model  for  Ang  shown  in  Figures  1-3  has  an  amplitude  error 
comparable  with  that  of  the  index,  about  30%.  However  the  leading  edges  of  the 
enhancements  are  located  quite  accurately;  the  error  is  about  6  hrs.  Details 
of  the  model  in  the  falling  portion  of  a  stream  are  less  reliable.  Here  the 
velocity  is  changing  slowly,  so  a  change  in  the  distribution  of  Ang  has  little 
effect. 

The  comparisons  of  Figures  1-3  are  summarized  in  Table  1.  Results  for 
both  the  midpoint  ("mid")  and  peak  ("pk")  velocities  are  shown.  Table  1  gives 
the  average  velocity  v,  the  standard  deviation  the  correlation  coefficient 
Pv  between  the  observed  and  model  velocities,  and  the  rms  difference  between 
these  velocities.  The  IPS  standard  deviations  and  the  p's  and  e's  are  all 
corrected  for  estimation  error  in  the  data.  For  comparison,  the  average  space¬ 
craft  velocity  was  499  km  s"^  and  its  standard  deviation  128  km  s  \  The 

correlation  coefficient  p  between  the  observed  and  model  scintillation  indices 

m 

is  also  shown.  It  was  calculated  after  removing  a  low  frequency  trend  from 
both  the  data  and  the  model  indices. 

Discussion 

Differences  of  detail  from  stream  to  stream  are  less  important  than  the 
general  character  of  the  Anp  distribution.  This  is  shown  in  Figure  4,  which 
gives  a  superposed  epoch  average  of  23  streams  aliqned  by  the  density  maximum 
at  the  leading  edge.  The  An&  profile  can  be  compared  with  epoch  averages  of 
other  plasma  parameters  or  used  as  a  prototype  weighting  function  when  inter¬ 
preting  IPS  data. 
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In  Figure  4  the  Ang  model  attains  a  maximum  twice  its  average  value  12  hours 

behind  the  density  peak.  It  decays  slowly  inside  the  stream,  falling  below  its 

average  value  after  four  days.  In  contrast,  the  proton  density  fluctuations, 

«n  ,  follow  the  bulk  density,  a  is  dominated  by  spatial  frequencies  ~ 10 km 
P  p 

Thus  the  difference  between  o  and  An  indicates  that  the  shape  of  the  fluc- 

np  e 

tuation  spectrum  changes  systematically  as  a  stream  passes.  At  the  density 
maximum,  the  spectrum  steepens  relative  to  its  pre-stream  configuration;  then  it 
flattens  inside  the  stream  before  relaxing  back  to  its  original  shape. 

In  Figure  5  we  divide  the  streams  into  two  groups  according  to  width.  The 
broad  streams  have  the  broader  Ang  enhancement,  as  though  high  speed  wind  were 
intrinsically  more  turbulent  than  low  speed  wind.  However,  our  data  from  high 
latitude  sources  do  not  support  this  interpretation.  These  sources,  which  are 
observed  through  faster  wind  than  in  the  ecliptic  (Coles  and  Rickett,  1976), 
show  scintillation  index  variations  with  elonaation  that  a^e  indistinguishable 
from  those  of  ecliptic  sources.  Furthermore,  there  is  no  apparent  difference 
between  index  data  of  1972  and  1974,  although  the  average  ecliptic  velocity 
differed  by  —  1 00  km  s (Feldman,  Asbridge,  Bame,  and  Gosling,  1978). 

The  extended  Ang  enhancement  in  the  broad  streams  appears  to  be  related  to 
the  internal  velocity  structure,  which  was  more  complex  than  in  narrow  streams. 

To  illustrate,  we  have  calculated  an  epoch  average  of  the  positive  velocity 
gradient,  VV+,  which  is  defined  to  be  zero  if  the  velocity  is  decreasing.  If 
a  constant  is  added  to  VV+,  as  in  Figure  5,  a  good  match  to  Ang  is  obtained. 
Apparently,  then,  An&  is  enhanced  in  regions  where  high  speed  wind  is  overtaking 
low  speed  wind.  The  enhancements  contribute  about  40%  of  the  total  micro¬ 
turbulence.  They  may  be  caused  by  a  process  unique  to  regions  of  increasing 
velocity,  or  they  may  simply  reflect  the  local  amplification  of  processes 
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occuring  everywhere.  If  indeed  the  strength  of  the  enhancements  is  determined 

by  the  velocity  gradient,  their  variation  with  solar  distance  may  be  slower 
_2 

than  r  ,  since  a  stream’s  velocity  profile  steepens  with  increasing  r.  Further¬ 
more,  the  enhancements  will  likely  be  less  pronounced  out  of  the  ecliptic,  where 
stream-stream  interactions  are  weaker. 

Figure  5  also  shows  the  proton  temperature,  T  ,  and  the  magnetic  field 
fluctuations,  <?g.  Both  Tp  and  Ang  are  enhanced  inside  the  stream,  but  their 
rising  and  falling  edges  do  not  coincide,  Og  and  Ang,  however,  appear  to  be 
closely  connected.  Og  is  measured  at  frequencies  ~ 10”**  km ,  but  the  distri¬ 
bution  of  field  fluctuations  at  higher  frequencies  is  probably  similar 
(Behannon,  1976).  If  so,  it  is  possible  that  Ang  is  physically  connected  to 
a  particular  type  of  field  activity,  for  example  fluctuations  due  to  hydromagnetic 
waves.  The  low  frequency  field  fluctuations  have  a  large  Alfvenic  component 
that  is  especially  pure  in  high  speed  wind  (Belcher  and  Davis,  1972).  However, 
Alfven  waves  cause  only  weak  density  fluctuations  and  thus  are  unlikely  to 
produce  appreciable  turbulence  directly  (Vlu  and  Huba,  1975).  On  the  other  hand, 
they  could  provide  a  driving  mechanism  for  the  production  of  high  frequency, 
turbulence-producing  waves  in  regions  of  overtaking  wind. 

Attempts  to  identify  such  high  frequency  waves  indicate  that  the  situation 
is  complex.  Unti  and  Russell  (1976),  for  example,  find  little  coherence  between 
high  frequency  plasma  flux  and  field  data,  indicating  that  relatively  pure 
waveforms  may  be  rare.  However,  Behannon  (1976)  tentatively  identifies  specific 
examples  of  ion  cyclotron  waves,  and  Neuaebauer,  Wu,  and  Huba  (1978)  show 
from  statistical  arguments  that  the  levels  of  the  flux  and  field  spectra  of 
Unti  and  Russell  (1976)  are  sometimes  consistent  with  the  predictions  for  magnet- 
osonic  waves.  But  the  latter  also  find  evidence  for  non-wave  field  and 


density  irregularities  which  are  maintained  by  internal  pressure  balance  and 
convected  rigidly  at  the  bulk  speed.  None  of  these  investigators  discusses 
regions  of  increasing  velocity  separately,  so  there  is  no  information  as  to 


the  characteristics  of  the  field  fluctuations  associated  with  the  An  enhancements. 

e 

We  plan  to  use  the  An  profile  of  Figure  4  as  a  known  weighting  factor  in 
an  attempt  to  determine  point  solar  wind  velocities  from  the  observed  sequence  of 
IPS  velocities  when  spacecraft  data  is  not  available.  We  will  concentrate  on 
non-ecliptic  regions  using  1973-75  data  from  high  latitude  sources.  This  investi¬ 
gation  should  also  reveal  the  extent  to  which  the  prototype  ecliptic  An  profile 
is  appropriate  for  high  latitudes. 
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Figure  Captions 


Figure  1:  Comparison  of  1973  IPS  midpoint  velocity  and  scintillation  index 
data  (the  connected  X's)  with  predictions  (the  solid  lines)  for 
Anfi  cc  and  the  final  Ang  model.  Linear  interpolation  in  the 

IMP  spacecraft  data  is  denoted  by  X's.  Units  for  velocity  are 

-1  -3 

km  s  ;  for  proton  number  density,  cm  ;  for  Ang  and  scintillation 

index,  arbitrary.  Vertical  lines  denote  Carrington  rotations. 

Figure  2:  Same  as  Figure  1,  but  for  1974. 

Figure  3:  Same  as  Figure  1,  but  for  1975. 

Figure  4:  Epoch  average  profiles  for  the  proton  velocity,  V;  the  proton 

density,  npi  the  microturbulence  model,  Ang;  and  the  low  frequency 
proton  density  fluctuations,  an  .  an  is  an  average  of  8  streams; 
the  others  are  averages  of  23  streams.  Each  parameter  (except  v) 
is  normalized  so  that  one  division  is  equal  to  the  parameter's 
mean  value,  which  is  denoted  by  a  bar.  Scales  have  offset  zeroes 
and  alternate  from  side  to  side. 

Figure  5:  Epoch  average  profiles  for  narrow  and  broad  streams,  showing  the 
proton  velocity,  v;  the  proton  density,  npi  the  microturbulence 
model,  Ang;  the  proton  temperature,  Tp;  the  low  frequency  magnetic 
field  fluctuations,  and  the  positive  velocity  gradient,  W+, 
to  which  a  constant,  c,  has  been  added.  The  last  three  are  over¬ 
plotted  with  Ang  (the  dashed  lines).  Each  parameter  is  normalized 
so  that  one  division  is  equal  to  the  parameter's  mean  value 
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(including  both  sets  of  streams),  which  is  indicated  by  a  bar. 
Scales  have  offset  zeroes  and  alternate  from  side  to  side. 
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"Solar  Cycle  Changes  in  the 
High  Latitude  Solar  Wind" 

S.  Ananthakrishnan,  W.A.  Coles,  J.J.  Kaufman,  and  B.J.  Rickett 


Abstract 

Measurements  of  the  solar  wind  speed  by  the  technique  of  interplanetary 
scintillation  are  presented  for  1971-79.  The  average  wind  speed  was  faster 
than  500  km/s  at  latitudes  above  +_  30°  for  most  of  1973-77.  This  fast  polar 
stream,  however,  was  observed  to  become  much  narrower  in  1978-79.  The 
narrowing  of  the  polar  streams  coincided  with  the  emergence  of  sunspots  at 
mid-latitudes,  with  the  start  of  the  new  solar  cycle,  and  with  a  corresponding 
contraction  of  the  polar  coronal  holes.  This  result  supports  the  idea  that  the 
solar  magnetic  field  controls  the  large  scale  structure  of  the  solar  wind. 
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It  has  long  been  known  that  the  solar  wind  cannot  be  spherically 
symmetric  although,  perforce,  this  simple  model  has  often  been  used.  The 
most  obvious  structures  are  the  fast  streams  which  are  now  known  to  have 
steep  spatial  gradients  both  in  longitude  and  latitude  (e.g.  Hundhausen 
1977).  We  present  here  measurements  of  the  latitude  structure  of  the  solar 
wind  speed  from  1971  to  1979,  made  using  interplanetary  scintillations  (IPS). 
The  results  show  that  persistent  high  speed  streams  emanate  from  the  poles 
of  the  sun;  and  that  these  streams  are  modulated  by  the  solar  cycle  exactly 
as  are  the  polar  coronal  holes. 

The  IPS  technique  has  been  discussed  by  many  authors  (Dennison  and 
Hewish  1967,  Armstrong  and  Coles  1972,  Watanabe,  Kakinuma,  Kojima,  and 
Shibasaki  1973).  Our  74  MHz  observatory  at  San  Diego  and  its  data  analysis 
are  described  by  Coles  and  Kaufman  (1968).  IPS  observations  in  the  ecliptic 
have  been  compared  with  direct  spacecraft  measurements  (Coles,  Harmon,  Sullivan 
and  Lazarus  1978).  That  work  gave  a  'calibration ‘  of  the  simple  interpretation 
which  we  adopt  here,  that  the  IPS  method  estimates  the  solar  wind  speed  at  the 
point  where  the  maximum  scattering  occurs.  We  make  IPS  observations  of  eight 
radio  sources  daily,  from  which  we  typically  obtain  useful  solar  wind  speed 
estimates  at  three  to  five  such  effective  locations.  The  heliographic 
latitude  of  these  locations  spans  60°S  to  70°N.  However,  the  coverage  is  not 
continuous;  northern  latitudes  (>  20°)  are  observed  from  March-July  and 
October-November,  while  southern  latitudes  (>  20°)  are  observed  June-July, 
each  year  (see  Fig.  1  of  Coles  and  Rickett  1979). 

The  solar  wind  speed  can  be  averaged  over  longitude  to  provide  the 
speed  versus  latitude  plots  of  Figure  1.  The  characteristic  l)-shape  has 
been  described  by  Coles  and  Rickett  (1976).  This  phenomenon  is  best 


described  as  evidence  for  persistent  fast  streams  from  the  north  and 
south  poles  of  the  sun.  These  polar  streams  are  the  solar  wind  response 
to  the  polar  coronal  holes,  in  just  the  same  way  that  ecliptic  streams 
originate  from  low  latitude  coronal  holes. 

Long-term  changes  in  the  polar  streams  can  be  seen  in  Figure  1. 

We  identify  three  phases. 

First,  1973-75  shows  stable  co-rotating  structures  with  two  "polar" 
streams  tilted  about  30°  from  the  rotation  axis,  as  revealed  in  Figure  2 
(Sime  and  Rickett  1976).  The  major  coronal  holes  at  this  time  can 
similarly  be  described  as  tilted  polar  holes;  Figure  3  shows  one  estimate 
of  the  large  coronal  hole  boundaries,  revealing  holes  well  aligned  with 
the  streams  shown  in  Figure  2.  Indeed  the  increased  average  speed  at 
the  equator  in  1974  (Figure  1)  is  due  to  these  streams,  as  discussed  in 
the  paper  by  Hundhausen  (see  this  volume). 

The  second  phase,  1976-77,  during  which  solar  activity  is  near 
minimum,  shows  polar  streams  that  have  widened  and  are  better  aligned 
with  the  rotation  axis.  This  gives  a  steep  average  latitude  gradient 
between  15°  and  30°  (over  10  km/s  per  degree)  comparable  with  instan¬ 
taneous  measurements  of  the  latitude  gradient  made  by  spacecraft  (Smith 
and  Rhodes  1974) . 

The  third  phase,  1978-79,  occurs  at  the  start  of  the  new  cycle  of 
solar  activity  and  shows  a  dramatic  decrease  in  the  speed  at  30-60°N, 
with  similar  changes  in  the  south.  In  other  words  the  polar  streams  are 
narrowing  at  the  same  time  that  the  new  sunspots  are  emerging  at  mid¬ 
latitudes.  This  change  evidently  occurred  at  the  same  phase  during 
the  previous  cycle.  The  change  that  we  observed  between  1977  and  1978 
was  seen  by  Hewish  and  Symonds  (1969)  between  1966  and  1967.  However, 
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they  observed  only  for  two  months  of  these  two  years  and  so  were  unable 
to  draw  any  strong  conclusions. 

The  presence  of  mid-latitude  sunspots  implies  that  bipolar  magnetic 
regions  are  dominating  these  latitudes,  restricting  the  area  of  open 
fields  (i.e.  the  polar  coronal  holes).  Evidence  for  the  reduced  area 
of  the  polar  holes  was  given  by  Sheeley  in  the  previous  paper  (see  this 
volume).  Figure  4  shows  a  comparison  of  average  wind  speed  at  0,  30°, 
60°N  with  the  area  of  the  north  polar  hole  from  1971-79.  Hole  area 
estimates  from  Kitt  Peak  He  10830  A  data  are  plotted  for  1974-78 
together  with  earlier  estimates  from  a  variety  of  techniques.  The 
slowing  of  the  wind  at  30°  and  60°N  in  late  1977  to  the  present  is 
closely  matched  by  the  shrinking  of  the  north  polar  hole. 

Figure  5  shows  the  comparison  as  a  plot  of  latitude  versus  time. 

The  continuous  line  in  the  north  and  horizontal  bars  in  the  south 
represent  the  boundaries  between  fast  and  slow  wind  from  six-month 
averages  of  wind  speed.  They  can  be  viewed  as  average  boundaries  of 
the  polar  streams  in  the  north  and  south.  The  dramatic  change  occurs 
between  the  period  1974  to  early  1977  and  the  period  late  1977  to  the 
present.  In  the  earlier  period  the  boundary  is  near  ±  30°,  showing  that 
the  fast  polar  streams  filled  50%  of  the  total  4-n  steradians;  while 
in  the  recent  period  the  boundary  moved  up  to  near  ±  65°,  showing 
that  the  polar  stream  now  fills  only  9%  of  4tt  steradians.  This  is 
similar  to  the  situation  in  1971,  although  our  data  coverage  for  that 
year  is  poor.  The  various  symbols  identify  the  polar  hole  boundaries 
(on  the  assumption  of  axial  symmetry).  Although  there  are  real  rotation 
to  rotation  changes,  the  long-term  behavior  can  be  described  as  a  60° 

22 


boundary  (13%  solid  angle  for  north  plus  south  hole)  in  1974-77, 
followed  by  a  retreat  to  70°  (6%  of  solid  angle)  in  mid- 1978 . 

We  may  define  a  solid  angle  expansion  ratio  between  the  (arbitrary) 

500  km/s  wind  speed  solid  angle  and  the  solid  angle  defined  by  the 

o 

He  10830  A  holes.  This  ratio  decreased  from  3.8  for  the  earlier  period 
to  1.5  for  1978.  Whereas  the  numerical  value  of  this  ratio  cannot 
easily  be  related  to  the  areal  expansion  factor  as  defined  for  the 
northern  polar  hole  by  Munro  and  Jackson  (1977),  the  change  in  our  ratio 
shows  that  the  areal  expansion  factor  of  the  polar  holes  decreased 
substantially  as  the  holes  shrunk. 

In  summary  we  have  observed  that  the  polar  streams  have  narrowed 
substantially,  at  the  same  time  that  the  polar  holes  have  narrowed. 

This  result  reconfirms  the  importance  of  the  solar  magnetic  field  in 
controlling  the  corona  and  solar  wind.  It  is  likely  that  this  changing 
three-dimensional  structure  determines  the  solar  cycle  modulation  of 
cosmic  rays.  From  the  viewpoint  of  SCADM,  it  is  clear  that  the  high 
latitude  solar  wind  shows  a  more  pronounced  solar  cycle  dependence  than 
does  the  ecliptic  solar  wind,  and  should  be  monitored  together  with  the 
sun  and  ecliptic  solar  wind.  Equally  important,  the  detailed  information 
on  the  latitude  structure  to  be  obtained  from  the  Solar  Polar  Mission 
must  be  viewed  in  the  conteA1.  of  the  long-term  evolution  as  presented 
here. 

We  thank  the  Atmospheric  Sciences  Section  of  the  N.S.F.  for  continued 
support  under  grant  ATM78-06770  and  previous  grants;  we  thank  AF6L  for  support 
under  contract  FI 9628- 77-C-0161 .  Wc  particularly  thank  our  colleagues  at  UCSD 
who  have  helped  with  the  project  over  the  years. 


Figure  Captions 

Figure  1.  Solar  wind  speed  measured  by  IPS  averaged  into  15°  intervals 
in  heliographic  latitude  for  the  years  1972-1978  and  January-June  1979.  The 
error  bars  are  +  twice  the  standard  deviation  in  the  average. 

Figure  2.  The  average  solar  wind  speed  in  15°  x  15°  bins  of  heliographic 
latitude  and  longitude  coded  by  the  grey  scale  as  indicated.  The  IPS  data 
were  averaged  for  Carrington  rotations  1612-1617  (March-July  1974). 

Figure  3.  Coronal  holes  indicated  by  a  contour  of  low  K-coronameter 

_  O 

brightness  (1.5  x  10”  pB ) ;  the  tick  marks  point  toward  low  brightness,  i.e. 
toward  the  holes.  Data  supplied  by  Drs.  R.  and  S.  Hansen  (H.A.O.). 

Figure  4.  North  polar  hole  area  as  percentage  of  solar  surface  (Sheeley, 
1979);  curve  is  drawn  by  eye.  Solar  wind  speed  from  IPS  observations  averaged 
over  six  months  and  15°  in  latitude  centered  on  the  latitudes  indicated. 

Figure  5.  Solid  curve  (in  the  north)  and  horizontal  bars  (in  the  south) 
indicate  the  contour  of  500  km/s  average  solar  wind  speed.  The  latitude 
boundaries  of  the  polar  holes  are  indicated  by  various  symbols:  4>  He  10830  ft, 
x  Bohlin  (1977),  Broussard  et  al .  (1978),  A  Fe  5303  A  (threshold  at  1.5  x 
10'6  of  disc  center  brightness,  data  from  Sacramento  Peak  Observatory). 
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